Dedicated to the memory of our friend Pierre Duclos 



o 

(N 



Conductivity and the Current- Current Correlation 

Measure 

Jean- Michel Combes, Francois Germinet, and Peter D. Hislop 

Abstract. We review various formulations of conductivity for one-particle 
Hamiltonians and relate them to the current-current correlation mea- 
sure. We prove that the current-current correlation measure for random 
Schrodinger operators has a density at coincident energies provided the en- 
ergy lies in a localization regime. The density vanishes at such energies and 

■ an upper bound on the rate of vanishing is computed. We also relate the 
current-current correlation measure to the localization length. 
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Q_|. 1. Current-Current Correlation Measure 

■ Higher-order correlation functions for random Schrodinger operators are 
essential for an understanding of the transport properties of the system. In 
this paper, we study the current-current correlation measure M(dE\,dE2) that 
describes the correlations between electron currents and is an essential ingredi- 
ent in the theory of DC conductivity for disordered systems. The connection 

' between the current-current correlation measure (referred to here as the ccc- 

measure) and the DC conductivity is expressed through the Kubo formula. 
We consider a one-particle random, ergodic Hamiltonian H w on £ 2 (Z d ) or 
I/"") ' on L 2 (R d ). In general H w has the form = ~(— iV - Aq) 2 + V per + V w , where 

QO . is a background vector potential, V per is a real-valued, periodic function, and 

V u is a real-valued, random potential. We always assume that H u is self-adjoint 
on a domain independent of u. We restrict H u to a box A C R d or Z d with 
Dirichlet boundary conditions. The finite-volume current- current correlation 
measures M^i, for a, f3 = l,...,d, are the point measures on M 2 defined as 
follows. Let E m (A) be the eigenvalues of H\ with normalized eigenfunctions 

4>m ■ Let V a H\ = i[H\,x a ] = p a — (A))a be the a th -component of the local 
velocity operator, where p a = —id/dx a is the momentum. The local velocity 
operator is independent of the randomness. We denote by v£J • the matrix 
elements of the a th -component of the velocity operator in the eigenstates (j)^ 
and <pf ] . That is, we define 

v£h = (^\v a H A 4% (1) 

We now define the finite-volume, current- current correlation measures for sub- 
sets A, cKby 

«*3(A I ,A,) = 1 ir E SAlr P> 

1 1 i,j:Bi(A)eAi;Bj(A)eA2 
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These set functions extend to point measures on 



written M„ A ' 
<?(Ai,A a ) 



Ma A d, are positive Radon measures on 



M, 



(A), 



The measures with a = j3, 
2 . We note that 

o Q (A2,Ai). We will denote by the positive measure 

( , L j.y± a ,a- The existence of the infinite-volume limits as A — > W d of these 
measures was proved in [25L I33L 134] . We denote the self-averaged, infinite- 
volume, positive, Borel ccc-measures by M a (dEi,dE 2 ), M a> p(dEi,dE 2 ), and 
M(dEi,dE 2 ). We remark that these measures can also be defined directly in 
the infinite volume using the trace-per-unit volume (see, for example, |25j). 

In order to relate the ccc-measures to the conductivity, we need to dis- 
cuss the Fermi distribution. The Fermi distribution n F (E;T) at energy E and 
temperature T > is defined for T > by 



n F {E-T)^(l + e^- E ^ T r\ 



whereas for T = 0, it is given by 



n F (E;0) = X(-oo,E F ](E). 



(3) 



(4) 



The density matrix n F (E F ; H u ) corresponding to the T = Fermi distribution 
is simply the spectral projector for and the half line (— 00, E F \. We will 
sometimes write P Ff for this projector. 

The Kubo formula for the AC conductivity at temperature T > and 
frequency v > is derived in linear response theory [8] (see also [3J, 128] ) by 
considering a time-dependent Hamiltonian with an external electric field at 
frequency v > 0. This Hamiltonian can be written as H u + £ ■ x cos ut, although 
it is more convenient to study this operator in another gauge. The conductivity 
relates the induced current to the electric field E(x,t) = £ ■ xcosvt. A formal 
calculation of the current to linear order in the electric field strength \£\ is 
sketched in section 14.11 The resulting Kubo formula for the AC conductivity 



is: 



lim lim 

e^O |A|^oo 



n F {E r ,T) - n F (E 2 ;T) 



Ei — E2 
npjE^T) - n F {E 2 -T) 

E\ — E 2 
n F (E + v;T)- n F (E;T) 



where 



8 e (s) 



5 e (E x -E 2 + v) M^(dE 1 ,dE 2 ) 
6(E 1 -E 2 + u) M a> p(dEi,dE 2 ) 
M a> p(dE,dE) . (5) 

(6) 



7T e 2 + s 2 



We mention that the ccc-measure has been computed in two models with the 
aim of computing the low frequency ac conductivity [161 I27j ; see also (281 
section 3.6]. 

We are interested here in the DC conductivity (T^F^ at zero temperature. 
Formally, this is obtained from ([5]) by taking the limits v — > and T — > 0. We 
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first take v — > in ([5]), to obtain 

<r a ,p(P> T) = - f dnF ^ T) M a ^(dE, dE). (7) 

Next, we take T — > 0. This results in a delta function —5{Ep — E). In order 
to evaluate the resulting integral, it is convenient to assume that the infinite- 
volume current-current correlation measure is absolutely continuous with re- 
spect to Lebesgue measure on the diagonal and that it has a locally bounded 
density m a> p(E,E). With this assumption, we obtain from ©, 

<J [ °p ) {E F )=m a ^(E F ,E F ). (8) 

Consequently, the diagonal behavior of the ccc-measure determines the DC 
conductivity. 

In this note, we investigate the existence of the densities m a ^ a {E,E) and 
upper bounds on the rate of vanishing for E in certain energy intervals. By ([8]), 
this implies the vanishing of the DC conductivity for energies in these intervals. 
The vanishing of the DC conductivity for energies in the localization regime has 
been proved by other methods, see section [H We also relate the localization 
length to the conductivity, and to the ccc-measure. 

We will consider the diagonal behavior of the ccc-measure at energies in 
the complete localization regime 3 . The complete localization regime is the 
energy regime discussed by Germinet and Klein [19J and is characterized by 
strong dynamical localization. We present precise definitions in the sections 
below. We will also discuss a possibly larger energy domain characterized by a 
finite localization length defined in (|45p . 

In order to formulate our results, recall that by the Lebesgue Differentiation 
Theorem, the positive measure M has a density at the diagonal point (E, E) if 
for any e > 0, with I e (E) = [E, E + e], we have that the limit 

e^o e 2 

exists and is finite. It is necessarily nonnegative. We say that the resulting 
density m(E, E) vanishes on the diagonal at the point (E, E) at a rate given 
by a function g > 0, with g(s = 0) = 0, if, in addition, for all e > small, we 

M(I e (E),UE)) 
< — v v y " < g(e). (10 
e z 

The boundedness of the current-current correlation measure on the diagonal 
means that the measure has no atoms on the diagonal. For comparison, the 
ccc-measure for the free Hamiltonian is easy to compute due to the fact that the 
velocity components V a H commute with the Hamiltonian. A simple calculation 
shows that 

M a>(3 (dE 1 ,dE 2 ) = n (E 1 )5(E 1 - E 2 )6 a ,pdE 1 dE 2 , (11) 

where uq{E) is the density of states of the free Laplacian at energy E. Con- 
sequently, the limit on the diagonal, as described in ([9]) does not exist and the 
DC conductivity is infinite at all energies. Boundedness may be obtained sim- 
ply for energies in a strong localization regime where ([22]) holds as stated in 
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Theorem [TJ If, in addition, a Wegner estimate (see [11 j . |12j ) holds, then the 
DC conductivity vanishes. 

We now state our main result on the vanishing of the ccc-density on the 
diagonal. 

Theorem 1. Let E E E . Then, the current- current correlation measure 
density m(E, E) exists and vanishes on the diagonal. If, in addition, the Wegner 
estimate (|39p holds on S , then for any bounded interval Iq C zf L with 
E G Jo, and for any < s < 1, there is a finite constant Cj 0tS < oo so that 

±M(I e (E),I e (E)) < C I(hS e|loge|s, (12) 

where I e (E) = [E,E + e}. 

We note that if, in addition to the Wegner estimate and the other hypotheses 
of Theorem [H we know that Minami's estimate is satisfied in a neighborhood of 
a given energy, then the rate of vanishing can be improved for certain intervals 
near the diagonal (see Corollary [1]) . We discuss this and some generalizations 
in section [2] after the proof of the main theorem. 

Remark 1. Suppose that Eq is a lower band edge at which the IDS exhibits 
a Lifshitz tail behavior. For example, if Eq is a band edge of the deterministic 
spectrum £ and E > Eq, the existence of Lifshitz tails [23, 29J means that for 
E sufficiently close to Eq, the IDS N(E) satisfies 



N(E) - N(Eq) < C Eo e ( B - B o) d/2 . (13) 

We write E(I e (EQ)) for the spectral projector for H u and the interval I t (EQ). 
It follows from (1130 that for e > small enough, we have 

E{Tv X QE(I € (E )) X Q}<C Eo e-^, (14) 

for some constant a > and where xo is the characteristic function for the 
unit cube about zero. We use this estimate in (j40j) below in place of the usual 
Wegner estimate (|39p . As a consequence, we obtain an exponential rate of 
vanishing on the diagonal. There is a finite, constant Do > so that for all 
e > small enough, 

\M(I e (E ),I e (E )) < D e- a ^ d/2 . (15) 

Of course, this only occurs at a countable set of energies (the lower band edges). 
This result, however, stating that, roughly, m{E + e, E) vanishes at an expo- 
nentially fast rate as e — > 0, is consistent with the Mott theory of conductivity. 

Remark 2. The optimal rate of vanishing in Theorem Q] is not known. We 
note, however, the relationship between the localization length (see section [3]) 
and the ccc-measure [51 section V]: 

e I Km,*®, (16) 

from which it follows that the ccc-density must vanish faster that 0(e) on the 
diagonal if the localization length is to be finite. We provide a proof of (j!6f) as 
part of Proposition [T] in section [3l 
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This leads us to our second result relating the localization length ^(A), for 
an interval A C R, to the ccc-measure. The localization length bounds the ccc- 
measure and the vanishing of the localization length implies the vanishing of 
the DC conductivity. For this result, we do not need the hypothesis of complete 
localization. 

Theorem 2. For E 6R, and for any e > small, define the interval I e {E) = 
[E, E+e]. Suppose that there is a constant < Me < oo so that the localization 
length 0,{I e {E)) < Me < oo, for all e > small. Then, the density of the ccc- 
measure m(E,E) exists on the diagonal. If, in addition, i{I t {E)) — > as e — > 0, 
then we have 

/.(g)) 

e->0 e 2 
Consequently, a^ DC \E) = 0. 

Finally, we mention the open question of the existence of a density 
m(Ei,E2) for E\ ^ E2 that we refer to as the off-diagonal case. This is im- 
portant in view of the formula for the localization length ()16p . and in order 
to control the density in a neighborhood of the diagonal. Unfortunately, the 
approach developed in this note does not seem to allow us to control the off- 
diagonal behavior. We mention that for lattice models on ^ 2 (Z rf ), Bellissard and 
Hislop [6] proved the existence of a density m a {E\, E2) at off-diagonal energies 
(Ei,E2) provided that 1) the energies lie in a region away from the diagonal 
E\ = E2 determined by the strength of the disorder, and 2) the density of the 
single-site probability measure extends to a function analytic in a strip around 
the real axis. 

2. Existence and Vanishing of the Density on the Diagonal in the 
Complete Localization Regime 

In this section, we give a simple proof of Theorem [1] on the existence of, 
and vanishing of, the ccc-density on the diagonal m(E,E) when E € E . We 
denote by C a generic, nonnegative, finite constant whose value may change 
from line to line. We work in the infinite-volume framework using the trace- 
per-unit volume and refer the reader to [8] and [5] for a complete discussion. 
Let (O, P) be a probability space with a Z rf -ergodic action r a : O -> fl, a 6 Z d . 
Let a G Z d — > U a be a unitary representation of Z, d on H. A covariant operator 
A on a separable Hilbert space H is a P-measurable, operator-valued function 
A = {A& I ui G ft} such that for a € lP ', we have {U a AU^ a ) w = A T _ aU . Let 
Xo be the characteristic function on the unit cube centered at the origin. The 
trace-per-unit volume of a covariant operator A, denoted by T(A), is given by 

T(A) = E{Tr( Xo ^ X o)}, (18) 

when it is finite. 

We recall that for Aj C M, with j = 1,2, the positive ccc-measure M a is 
given by 

M a (Ai,A 2 ) = T(V a HE(A 1 )W a HE(A 2 )) 

= lim JL Tr X AV Q ^(A 1 )V a ^(A 2 ) X A. (19) 

|A|->-oo |A| 
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We recall from [8j section 3.2] that K 2 is the space of all measurable, covariant 
operators so that 

E{||A X o||l} < 00. (20) 
We refer to [8 for further definitions and details. Because the spectral projector 
E(A) G K, 2 , the cyclicity of T permits us to write 

T(V a HE(A 1 )V a HE(A 2 )) = T(E(A 2 )V a HE(A 1 )V a HE(A 2 )). (21) 

Let Xx be a bounded function of compact support in a neighborhood of 
x G M d . Let || ■ ||p denote the p^-trace norm for p > 1. We give a characterization 
of the region zF L of complete localization given in [18j and |19j . 

Definition 1. An energy E G £ belongs to S if there exists a open neighbor- 
hood Ie of E such that the following bound holds. Let Ie denote all functions 
f G with supp/ C Ie and sup tgR |/(*)| < 1- Then, for any < s < 1, 

there is a finite constant C s > so that the Hilbert- Schmidt bound holds 

sup H\\Xxf{H u )xy\\l} < C s e^ x -^ s , for all x,y G R d . (22) 
fei E 

Remark 3. This estimate can be improved to s = 1 by the method of fractional 
moments [21 14] . 

Proof of Theorem [TJ 1. We begin with formula (|2T|) and take Ai = A 2 = 
I t {E) = [E,E + e]. We write the velocity operator as a commutator 

V a H = i[H,x a ] = i[H - E,x a }. (23) 

For brevity, we write P e = E(I t (E)). We then have 

M a (I e (E),I € (E)) = T(P e V a HP e V a HP e ) 

= -T(P e [H-E,x a ]P,[H-E,x a ]P e ). (24) 

Next, we expand the commutators and obtain four terms involving P t . This 
is facilitated by the introduction of a a function / e (s) = (s — E)xi t (E){ s ) / e -> so 
that |/ e | ^ 1. We obtain a factor of e from each projector so that 

M a (Ie(E),I e (E)) = -e 2 {I + II + III + IV), (25) 

where we have: 

/ = T(f t (H)x a f t (H)x a P t ); (26) 

// = -T(f e (H)x a P e x a f e (H))- (27) 

/// = -T(P e x a f?(H)x a P e ); (28) 

IV = T{P e x a f e {H)x a f e {H)). (29) 

2. In order to estimate these four terms for E G zF L , we need the following 
bounds. Let Fa^H^) denote the operator x a fe{Huj)Xa- The bound in Defini- 
tion Q] implies the following estimate for any integer q > and s as in definition 

m 

supEiWiF^H^fiHMll} < C(s,q,d)T((q + d)/s) 2 , (30) 
felt 

where T(t) is the gamma function and C(s, q, d) depends on s through the con- 
stants C s and a s as in (|22"j) and grows linearly in q. To prove this bound, we 
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use a partition of unity on W 1 given by Y^k& d Xfc = 1 an d we define opera- 
tors A ki = Xkf{H u )xe. and B M = Xk{x a f{H^)x a ) q xi- We then estimate the 
Hilbert-Schmidt norm in (|30p by 

UF^H^fiH^xoWl = \\xof(H u ){F a , e (H u ))*f(H u ) X o\\i (31) 

< II^Omlbll^olbll-Bm^ll- (32) 



fit 



kez d 



Since ||-B TO ^[| < C(g)||m||' ? ||^||' 3 , where C(q) is linear in q for integer q, we obtain 
from (|22p and (|31|) . and the integral representation of the gamma function 

E{||(F a , e (^))"/ 2 /(^)xo||i} < C( S ) ( £ |m|V-*IHI' J (33) 

< C(*,g,d)r((g + d)/ a ) 2 , (34) 
proving the estimate ([30]) , We will also use the operator bounds 

\\f(H)\\<l, |/|<1, (35) 
and by a simple trace class estimate, 

T(P e ) < C, (36) 

for some finite constant C > 0. 

3. We now bound the term / in (I26p as follows, recalling that \f t \ < 1 and that 
Pefe{H u ) = f e {H w ): 

\I\ < \E{Tr X ofe(H)x a f e (H)x a P eX o}\ 

< T(P e ) 1/2 (E{\\F a ,e(H.)fe(HM\ 2 2}) 1/2 

< r(Pe) V2+l/4 (E{||F a , e (^ 2 / e (^) Xo |!i}) 1/4 (37) 



Bounds (|30j) . together with ([35^ - (|36]) . show that the term in (j37|) is uniformly 
bounded as e — > 0. It is clear that terms 11 — IV are bounded in a similar 
manner. As a consequence, it follows from these bounds and (|25j) that 



h wa< c<m (38) 

e^O e 2 

This implies that the density m(E,E) exists for E £ E . 

4. The rate of vanishing can be calculated with more careful estimates. We use 
the Wegner estimate in place of (|36h in order to obtain more powers of e. The 
Wegner estimate for infinite-volume operators has the form 

T{E{J)) = E{Tt X oE(J)xo} < C \J\, (39) 



for any subset Jcl Turning to term I in (137p . we iterate the argument n 
times, recalling that Pef^H^) = f^H^)-, and obtain 

\I\ < r(i> )*+•••+£ (EmF^H^) 2 "- 1 f e (H u )xo\\ 2 2 })^ ■ (40) 
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We use the strong localization bound (|30p in (|40p . together with the Wegner 
estimate (|39]) for I e (E), and Stirling's formula for the gamma function, to obtain 



|/| < C(s,d)e 1 -^ (2 n )i (41) 

We choose n so that 2™ ~ | log e| producing the upper bound e|loge|s. The 

remaining terms can be estimated in a similar manner. □ 

Remark 4. We note that if we consider intervals of the form 

I e = [E,E + e] x [E-e,E] CR 2 , (42) 



we can improve the rate of vanishing. This is a consequence of the analysis in 
[28J in the context of Mott's formula for the AC conductivity. It requires the 
Minami estimate in addition to the Wegner estimate. We recall that Minami's 
estimate for local operators is a second-order correlation estimate. We 
assume that the single-site probability measure has a bounded density p. Let 
E\(A) denote the spectral projector for and the interval A. The Minami 
estimate is: 

E{Tr£ A (A)(Tr£ A (A) - 1)} < (\\p\UAl\A\) 2 . (43) 

This was proved by Minami |31j for lattice models on and recent, simplified 
proofs have appeared in [7J [9l 124] . More recently, the Minami estimate was 
proved in [10J for certain Anderson models in the continuum for an interval of 
energy near the bottom of the deterministic spectrum. 

Corollary 1. For any energy E E S for which both the Wegner and Minami 
estimates hold in an interval Iq containing E, we have 

L M ([E,E + e],[E-e,E]) <C Io e 2 |loge| d+2 . (44) 

Proof. As mentioned, this result is a consequence of the analysis in |28j of 
the Mott formula. We go back to the decomposition given in (j25j) of that paper, 
and estimate each term using [28, Theorem 4.1]. Applying }28[ Theorem 4.1] 
directly to bound M([E, E + e], [E - e, Ej), we obtain (jSJ). □ 

3. Relation to localization length 

The second moment of the position operator is used in the following co- 
variant definition of the localization length due to Bellissard, van Elst, and 
Schulz-Baldes [5]. 

Definition 2. The localization length for an interval A C R is 

£(A) 2 =hm i f T(E(A)\x(t)-x\ 2 E(A)) dt, (45) 

where x(t) = e~ ltH xe ttH . 

If A C R is contained in the complete localization region it is easily seen 
from (|22p that the localization length is finite but one expects that the finite- 
ness of the localization length holds in a much larger energy domain. Notice, 
however, that as shown in [5], the finiteness of ^(A) implies that the spectrum 
of H u in A is pure point. This does not a priori imply estimates such as (|22l . 
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The ccc-measure on the diagonal and the localization length are very closely 
related. Recall that the velocity operator is = i[H,x]. 

Proposition 1. Let E(dX) be the spectral family for H^. The ccc-measure is 
related to the localization length by 

£{A) 2 = 2 [ [ T(E(dfi)VHE(d\)VHE(dfj.)) _ 
■I ./a (A-//) 2 
Furthermore, we have the bound 

T{E{A)VHE(A)VHE{A)) 



< 1(A) 2 . (47) 



Proof. We use the Fundamental Theorem of Calculus to write 

r-t 

-isH -t r AsH 




x(t)-x = - e- tsH Ve lsH ds. (48) 
Jo 

Taking the absolute square, we obtain 

\ x (t) - x\ 2 = f dw f ds e~ isH VHe^ s " w)H VHe iwH . (49) 
Jo Jo 

Using the spectral family for H in the form J K E(dX) = 1, we obtain 
T(E(A)\x(t) - x\ 2 E(A)) 

dse~ is ^- X) jf dwe iw(p - x) T(E(dfi)VHE(d\)X7HE(dfi))^ . 

(50) 

We first perform the integration over s and w. We next take the time average. 
Finally, taking the limit T — > oo, we obtain (|46p . 

In order to obtain (|47|) . we write the projectors in the trace-per-unit volume 
on the left of (0ZJ) as E(A) = J A E(dfx) and note that (A - n) 2 < \A\ 2 for 
X,fj,eA. □ 

The simple proof of Theorem [2] follows from (|46p . For Bel, and e > 
small, we have 

%M(I e {E)J e {E)) = % I I T(E{dn)VHE{d\)VHE{dn)) 
e e J Ji e (E)xic(E) 

< t(h(E)) 2 . (51) 

Taking E\ = E2, we obtain the proof of the second part of the theorem. 

We next present an explicit formula relating the localization length £(A) to 
the second moment of the position operator. 

Proposition 2. We have the following identity: 

\i{A) 2 = T(E(A)\x\ 2 E(A)) - [ T(E{dX)xE(dX)xE(dX)) > 0. (52) 
2 J a 

where E(dX) is the spectral family for H^. Consequently, if the second moment 
of the position operator T(E(A)\x\ 2 E(A)) on the right side of \5S\) is finite, 
then the localization length is finite. 
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Proof. The proof of the equality in (|52p follows from (|46p by expanding 
the commutator in the definition of VH. The nonnegativity is obvious and this 
proves the second part of the proposition. □ 

We remark that the formulas in (|45p and in (|52|) are manifestly covariant 
with respect to lattice translations. It is known [5] that £(A) < oo implies 
that the spectrum of H w in A is pure point almost surely. The finiteness of 
the term T(E(A)\x\ 2 E(A)) should be compared with the finiteness condition 
of the Fermi projector given in (|57p for the continuum model and in (|65p for 
the lattice model. 



4. Conductivity for Random Schrcidinger Operators 

In this section, we discuss the various notions of conductivity that have 
occurred in recent literature and provide a justification for the calculation of 
the conductivity presented in section [TJ As a starting point, we begin with the 
adiabatic approach to the Kubo formula as presented in [§]. We then consider 
related results of Frohlich and Spencer (1983) [TF], Kunz (1987) [30], Bellissard, 
van Elst , and Schulz-Baldes (1994) [5], Aizenman-Graf (1998) [3], and Nakano 
(2002) [32j. We write v a or V a H for the a th -component of the velocity operator 
v a = V a H = i[H,x a \. 

We mention another derivation of the Kubo-Streda formula for the trans- 
verse conductivity in two-dimensions for perturbations of the Landau Hamil- 
tonian at zero frequency and temperature is presented in [131 section IV]. This 
derivation does not require the use of an electric field that is adiabatically 
switched on. On the other hand, the Fermi level has to be restricted to a 
spectral gap of the Hamiltonian, rather than to the region of localization. 

4.1. The Adiabatic Definition of the Conductivity. 

4.1.1. Kubo formula through linear response theory. In the paper of Bouclet, 
Germinet, Klein and Schenker [8], the total charge transport is calculated in 
linear response theory in which the electric field is adiabatically switched on. 
This provides a rigorous derivation of the DC Kubo formula for the conduc- 
tivity tensor. A noncommutative integration approach is presented in |15L I14j 
which somewhat simplifies manipulation of operators. The random, Z rf -ergodic 
Hamiltonian H w on L 2 (R rf ) has the form 

H u = (iV + A u) f + V ul , (53) 

where (Aj ; KJ are random variables so that this operator is self-adjoint. We 
add a time-dependent, homogeneous, electric field £{t) that is adiabatically 
switched on from t = —oo until t = 0, when it obtains full strength. Such a 
field is represented by £{t) = e 71 £, where £ is a constant and t- = min (0, t). 
The usual gauge choice is the time-dependent Stark Hamiltonian given by 

H u (t) = H u + £{t) ■ x. (54) 

This Hamiltonian is not bounded from below. As is well-known, and used in 
[8] , one can make another choice of gauge to eliminate this technical problem. 
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Let F(t) be the function given by 

F(t) = f £(s) ds, (55) 

J — oo 

and note that the integral converges provided rj > 0. 

We now consider another Hamiltonian obtained from in (|53p by a time- 
dependent gauge transformation using the operator G(t) = e iF ^' x : 

H u {t) = G(t)H u G*(t) = (iV + K + F(t)) 2 + V u . (56) 

This operator is manifestly bounded from below almost surely provided is 
also. Furthermore, the two Hamiltonians H u (t) and H u (t), given in (p3|) . are 
physically equivalent in that they generate the same dynamics. If tj)t solves the 
Schrodinger equation generated by H w (t), then the gauge transformed wave 
function G(t)*i/jt solves the Schrodinger equation associated with H u (t). Con- 
sequently, we will work with H w (t). 

In order to describe the current at time t = 0, the system is prepared in an 
initial equilibrium state £ w at time t = — oo. The initial state is usually assumed 
to be the density matrix corresponding to the Fermi distribution at temperature 
T > 0. If T > 0, the state with Fermi energy Ep is given by £ w = npiH^; T), 
where np(E;T) is given in If T = 0, the initial equilibrium state is the 
Fermi projector £ w = np(H u} ;0) = X(-ca,E F ](Hu))- We will often write Pe f for 
this projection. More general states are allowed, see [U section 5]. The crucial 
assumption for the T = case is 

E{|||:e|P£ fXo ||2}<oo, (57) 

which holds true whenever the Fermi level lies in an interval of complete local- 
ization zF L [20J. If the initial state is given by a Fermi distribution, together 
with (|57h if T = 0, the authors prove the existence of a unique time-dependent 
density matrix Pu(t) solving the following Cauchy problem: 

id t Puj{t) = [H u (t),p u (t)]* 
lim p u (t) = Ce f , (58) 

t— y — oo 

in suitable noncommutative L 1 and L 2 spaces (see |14L 115] for the construction 
and use of the more natural IP spaces defined over the reference von Neumann 
algebra of the problem, and the corresponding Sobolev spaces). In particu- 
lar, the commutator in (158j) requires some care in its definition when is an 
unbounded operator. The subscript star reminds the reader of these noncom- 
mutative integration spaces. 

In [8], the authors also prove that in the T = case, the density matrix 
Pw(t) is an orthogonal projection for all time. Given this density matrix p w (t), 
the current at time t = is defined to be: 

J( V ,£;E f ) = T(v(0)pU0)), (59) 

where the modified velocity operator is v(0) = i[H u , x] — 2F(0) = v — 2F(0). At 
time t = — oo, the system is assumed to be in equilibrium so the initial current 
is zero. The current J(rj,£; Ep) is the net current at time zero obtained by 
adiabatically switching on the electric field. 
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Let U(s,t) be the unitary propagator for the time-dependent Hamiltonian 
(|56|) . We denote by Ce f (t) the time-evolved operator obtained from the Fermi 
distribution rip (_ff w (r) ; T) , as in ([3]) for T > 0, or as in (j3J), for T = 0. The 
components of the current are given explicitly by 

J Q (V,£;E F ) = -T\ f dre^v a (0) (U(0,r)(i[£ ■ x,£ Ef (t)})U(t,0))^ 
I J — oo 

(60) 

Linear response theory now states that the conductivity is the constant of 
proportionality between the current J and the electric field £ , that is, that 
Ohm's Law holds: J = a£. This is a tensorial relationship. Formally, one 
should expand J(r],£]Ep), given in (j60l) . about £ = 0. The components a a p 
are obtained via the derivative: 

/ dJ a (n,£;E F ) \ 
a a p{-q,E F )=y 7^- J ■ (61) 

The question of the limit £ — > is addressed in the following theorem. 

Theorem 3. [8j Theorem 5.9] For rj > 0, the current defined in \59\) is dif- 
ferentiable with respect to £j at £ = 0, and the conductivity tensor is given 
by 

a af3 (rr,E F ) = -Ty° dr v a (e"^^.^])^) \ . (62) 

Once again, the subscript star in (|62p indicates that operators and product of 
operators are considered within suitable noncommutative integration spaces. 
Formula (|62p is valid for a large family of equilibrium initial states including 
the Fermi distributions. The analogue of [5j Eq. (41)] and [351 Theorem 1] then 
holds. 

Corollary 2. pS Corollary 5.10] Assume that £(t) = Ree iut , then the 

conductivity o~jk(ri;^;u) at frequency v is given by 

a jk (n;(;u;0) = -T {v a {id + r, + iv)' 1 (fyC)}} » (63) 

where L\ is the Liouvillian generating (|58p . and dpQ = i[xp,^E F ]- 

For T = 0, following ideas of [5], the authors recover the Kubo-Stfeda 
formula for the conductivity tensor. The distribution £ Ef is given by the Fermi 
projector as in @. In this case, we can take the limit rj — > and obtain an 
expression for the DC conductivity. 



Corollary 3. [8j Theorem 5.11] Under the hypothesis (51) on the Fermi pro- 
jector P Ff , the conductivity tensor is given by 

a^(E F ) = a a p = -%T {Pe f [[xoc,Pe f },[xp,PeM ■ ( 64 ) 

The tensor is antisymmetric, so that a aa = 0, for a = 1, . . . , d. If, in addition 
to (57\ ), the magnetic vector potential in |53|) vanishes, A u = 0, so that the 
Hamiltonian is time-reversal invariant, all components of the conductivity 
tensor vanish. 
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Hence, all components of the DC conductivity vanish for time-reversal invariant 
systems (Hamiltonians (|53p with ^4 = 0) provided the Fermi projector satisfies 
the decay hypothesis (|57p . We note that for lattice models, this condition is 

m Yl n 2 \(°\ p ^\ x )\ 2 \ <00 < ( 65 ) 

where \y) denoted the discrete delta function at y € Z d . 

For systems with nontrivial magnetic fields, the off-diagonal terms may be 
nonzero. Of special interest is the two-dimensional case with a nonzero, constant 
magnetic field. In this case, the integer quantum Hall effect is the fact that the 
off-diagonal term <Ti2 is an integer multiple of e/h when the Fermi energy Ep 
is in the region of localized states (see [5J, [3] for lattice models and )21[ [22 J for 
models in the continuum). 

4.1.2. Relation to the ccc-measure. The formula derived from rigorous linear 
response theory (j62|) can be manipulated in order to derive the T = formula 
©. We write ([62]) as 

/•oo 

a a>p (ri;E F ) = -2i / dr e^ T T(v a e- irH [xp, Pp F ]e lrH ). (66) 
Jo 

We use the spectral theorem H = J XdE\, and the identity 

[xp,f(H)]=i[ f(s) f du e-^ H vpe- luH . (67) 

Jr JO 

Recalling that the Fermi projector Pp F = njr(i?;0), we formally write the 
trace-per-unit volume as 

T(v a e- H [x,,P EF ]e- H ) = J J ^^^^T(v a dE u v,dE x ) (68) 

Substituting (|68p into (|66p . and performing the time integration with use of 
definition ([6]), we obtain 

a a AE F ) = lim/ I MX; 0) ~ Hf{v] 0) 5 e (X - v)T(v a dE u v p dE x ) 

= lim / / "HA;0)-nH*;0) 5e{x _ v)M {dUjdX)] (69) 

J \ — v 

where M a ^ is the ccc-measure. 

4.2. Results on DC Conductivity. We review several papers concern- 
ing the conductivity in one-particle systems. Typically, these authors assume 
a formula for the conductivity and then show that under a condition such as 
([57]) or ([65]) the DC conductivity vanishes. 

4.2.1. Frohlich- Spencer and Decay of the Green's Function. In their funda- 
mental paper of 1983, Frohlich and Spencer |17j proved the absence of diffusion 
and the vanishing of the DC conductivity for the lattice Anderson model. The 
authors assumed a form of the Kubo formula that expresses the DC conduc- 
tivity of a gas of noninteracting electrons at T = and Fermi energy Ep in 
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terms of the Green's function. Let p(E) > be the density of states for H u . 
For a, f3 = 1, . . . ,d, the conductivity tensor is given by 

vj$(E F )p(E F ) = lim ^ V x a xpE{\G(0,x;E F +ir})\ 2 }. (70) 
" J7-»0 air i — ' 

x&L d 

Frohlich and Spencer showed that if Ep lies in an energy interval for which the 
multiscale analysis (MSA) holds, then cr^(Ep) = 0, assuming p{Ep) > (see 
[261136]). 

Theorem 4. [17] Suppose the Hamiltonian = L + V w on £ 2 (7, d ), where L 
is the lattice Laplacian (in particular, the magnetic field is zero). Suppose that 
bound |7i| ) holds for Lebesgue almost all energy in (a,b). Then if Ep £ (a,b), 
the DC conductivity defined in J70|) at Fermi energy Ep vanishes. 



The key MSA bound on the Green's function is that, with probability one, 
Bwp\G u (x,y;E + iri)\ < de'" 1 '^. (71) 

r]>0 

The exponential decay of the Green's function then implies that c^fa = 0. 
Formula (|70p is equivalent to (|69p . The following derivation is presented, for 
example, in [3]: 

2 

^- V x a xpE{\G(0,x;E F +ir))\ 2 } 

IT i — ' 

2 

= -1-E{(0\[x a , {H-E F + ir?)" 1 ]^, (H-Ep- ir,)~ l ]\Q)} 

IT 

2 

= —T{R(E F - irj)v a R(Ep - irj)R(Ep + iri)vsR(Ep + irj)} 

IT 

= J J 8 v {X-Ep)5 r) {p-E F )M 0l ,p{dX,dv). (72) 

If, as in section [I] we assume the existence of a ccc-density m a g, we may 
compute the limit rj — > in ([72]) . We then obtain a^(Ep) = m a ^(Ep, Ep). 
So, by ([8]), this gives the same conductivity as defined in ©-(UJ). 

4.2.2. Kunz. Kunz [30] considered the two-dimensional Landau Hamilton- 
ian with a random potential on L 2 (R 2 ) so that H w = (1/2) (iV + A ) 2 + V u , 
where Aq = B(x2,0) and the velocity operators Vj = i[H u ,Xj]. He assumed 
that the conductivity is given by a Kubo formula that he wrote as 

1 f'°° 

ai, 2 (E F ) = lim - / e- et T(PE P lvi,e ttH »v 2 e- itH «]P Ep ) dt, (73) 
e^O ie J 

and for the diagonal terms 

1 f'°° 1 
vjjiEp) = lim - / e-^T{Pp F [v^e^Vje-^PEp) dt + -T(P Ef ), (74) 

for j = 1,2. Note that ([73]) and ([71]) follow directly from the above Kubo 
formula ([62]) . with £e f = P Ep , after integrating by part (and a change of 
variable t — > —t). Assuming band-edge localization for H u , he shows that if 
the disorder is weak enough relative to B so that there is a gap between the 
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Landau bands and if the Fermi energy Ep lies in the gap between the n and 
(n+l) s *-bands, then the transverse conductivity ai^Ep) is a universal multiple 
of (n + 1). He also provides arguments for the localization length to diverge in 
each Landau band. 

4.2.3. Bellissard, van Elst, and S chultz-B aides 's Work on Lattice Models. 
Bellissard, van Elst, and Schulz-Baldes [5] derived (|64p in a one-electron model 
using a relaxation time approximation. The one-particle Hamiltonian differs 
from (f54"|) in that a time-dependent perturbation W co u(t) is added that mimics 
a dissipation process. The interaction has the form W co u(t) = J2kez Wk8(t— t^). 
The ordered collision times tk are Poisson distributed so that = t& — t^-i 
are independent, identically distributed random variables with an exponential 
distribution and mean collision time r = E(rfc). The amplitudes are the 
collision operators that are assumed to commute with H and be random oper- 
ators. The model is discrete. They computed the time-averaged current using 
this evolution and found 

2 i 

•WW = Yl £iT([xi,n F (H; T)] 6 + f _ CH _ g ^ [x,H]). (75) 

i=l 

where jCh(A) = i[H,A] is the (bounded) Liouvillian. The authors then ne- 
glected the 8 ■ V term in the resolvent appearing in (|75p and took the limit 
5 —> 0. This exists provided the collision factor f > 0. This factor is pro- 
portional to 1/r, the relaxation time. Upon differentiating with respect to the 
electric field, they obtained 

o- a ./3 = T([x a ,n F (H; T)] _ [xp, H}). (76) 

r - Lh 

In certain cases, the temperature T can be taken to zero and the relaxation 
time can be taken to infinity so f — > 0. For example, they proved that for 
the Landau Hamiltonian with a random potential in two-dimensions, the off- 
diagonal conductivity crf^ (E) agrees with Kubo-Stfeda formula (i64"|) and is 
a constant on energy intervals where (|65|) holds. Moreover, it is an integer 
multiple of e 2 /h. 

4.2.4. Aizenman-Graf's Work on Lattice Models. Stimulated by the inte- 
ger quantum Hall effect, Aizenman and Graf [3] considered the analog of the 
randomly perturbed Landau Hamiltonian on the lattice. In this case, the DC 
conductivity for the system with a constant magnetic field is given by the Kubo- 
Stfeda formula: 

o-°$(E) = iTr{P E [[x at P E ], [x0,Pe]}}, (77) 

provided the operators are trace class. The exponential decay of the kernel of 
the Fermi projector (|84p implies (|65p . 

4.2.5. Nakano's Result for Lattice Models. The vanishing of the diagonal 
terms of the DC conductivity tensor for random, ergodic Schrodinger operators 
on the lattice Z, d was also proven by Nakano [32J provided the Fermi projector 
satisfies ([65]) . 

Nakano defined two (scalar) conductivities associated with the xi-direction 
(any coordinate direction can be used) differing in the order in which certain 
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limits are taken. Let H w £ = H u + Sx±, corresponding to the Anderson Hamil- 
tonian H w with an electric field in the xi-direction. The current of the system 
at zero temperature and at time t is given by 



Ji(t;S;E F ) = T(e UH ^ Vl e- UH ^P EF ), (78) 
where v\ = i[H w ,xi\. The corresponding conductivity is given by 

a a (t; Ep) = lim -J x {t; £; E F ). (79) 

consistent with linear response theory. The first definition of the conductivity 
a a (Ep) is the time average of the conductivity a a (t;E F ) defined in ([79]) : 

a a (E F ) = lim i f dtax(t;E F ). (80) 

The second definition of Nakano replaces the time-average by an Abel limit: 

1 f°° 

a h (E F ) = lim lim - / dt5e~ St Jx(t]£]E F ). (81) 

This is the same as the result derived in [5] using the relaxation time approxi- 
mation. 

Theorem 5. Under the assumption \65}) on the Fermi projector Pe f , with \x\ 
replaced by \x\\, we have 

a a {E F ) = a b {E F ) = 0. (82) 



4.3. Localization and the Decay of the Fermi Projector. When does 
the hypothesis (|57"|) for continuum models, or (|65j) for lattice models, hold? 
Aizenman and Graf [3] used the method of fractional moments, developed in 
|4j and [T] for lattice models on £ 2 (Z d ), to prove the decay of the Fermi projector 
(|65p provided the Fermi energy is in the strong localization regime. 

The basic hypothesis is that the fractional moment estimate holds: 

supK{\G(x,y;E + in)\ s } <C s e~ s ^ x ~ yl , x,y£Z d , (83) 

7)>0 

uniformly for n > 0, for some < s < 1 and \i > 0. If (|83j) holds at energy E F , 
then they prove that the kernel of the Fermi distribution at T = satisfies 

E{<s|P^|y)}<C7 e-^l. (84) 

As pointed out by Aizenman and Graf, this exponential decay estimate on the 
Fermi projector requires only that E F lie in a regime of energies for which the 
fractional moment bound (|83h holds. It does not require that (I83h hold at all 
energies below E F . The fractional moment bound (j83|) is known to hold in the 
strong localization regime or at extreme energies, cf. pQ, [4]- 

At last, we mention that fast decay of the kernel of the Fermi projection at 
a given energy E turns out to be equivalent to E G zf L |20j. 
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